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putational simulations and molecular line observations of the Taurus molecular 



I cloud. A new analysis method is presented to evaluate the degree and angle of 

I velocity anisotropy using spectroscopic imaging data of interstellar clouds. The 

efficacy of this method is demonstrated on model observations derived from three 
dimensional velocity and density fields from the set of numerical MHD simula- 
r-| I tions that span a range of magnetic field strengths. The analysis is applied to 

^1 ^^CO J=l-0 imaging of a sub- field within the Taurus molecular cloud. Velocity 

O I anisotropy is identified that is aligned within ~10° of the mean local magnetic 

^ I field direction derived from optical polarization measurements. Estimated val- 

ues of the field strength based on velocity anisotropy are consistent with results 
from other methods. When combined with new column density measurements for 
^ I Taurus, our magnetic field strength estimate indicates that the envelope of the 

^ I cloud is magnetically subcritical. These observations favor strong MHD turbu- 

O I lence within the low density, sub-critical, molecular gas substrate of the Taurus 

' cloud. 
O , 

OO , Subject headings: ISM: clouds - ISM: magnetic fields - ISM: kinematics and 



dynamics - ISM: individual (Taurus Molecular Cloud) - physical data and pro- 
cesses: MHD; methods: data analysis 



1. Introduction 

Dense, interstellar molecular clouds offer a unique and valuable laboratory to investigate 
magneto-turbulent phenomena. These clouds are expected to be fully turbulent systems with 
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a very large dynamic range between driving and dissipation scales. The degree of magnetic 
coupling to the turbulent flows has important implications for the nature of gas dynamics 
and star formation within molecular clouds. A strong, well coupled field can affect the star 
formation efficiency in a cloud by reducing the amount of material that is susceptible to 
gravitational collapse and star formation, and also affect the scale at which collapse occurs 
(Mouschovias 1976; Vazquez-Semadeni etal 2005). Magnetic fields also strongly affect the 
degree of gas density compression in shocks. Such shock-generated density perturbations 
may provide the seeds of protostellar cores and protoclusters. Given the potential impact of 
the magnetic field on the gas dynamics of molecular clouds, it is imperative to measure (or 
estimate) magnetic field strengths and to develop accurate descriptions of magnetohydrody- 
namic (MHD) turbulence under conditions applicable in star-forming clouds. 

Goldreich and Sridhar (1995, hereafter GS95) developed a theory for strong, incompress- 
ible, MHD turbulence that provides definitive predictions of the spectrum and anisotropy 
of velocity fields. Wave-wave interactions are expected to shear the Alfven wave packet in 
the plane perpendicular to the mean field. Correspondingly, wave energy is more efficiently 
redistributed to smaller scales in the direction perpendicular to the field than through the 
cascade parallel to the field. GS95 propose that a critical balance is achieved between non- 
linear interactions and wave propagation, such that the time scales to transfer energy along 
the two directions are comparable. 



where X\\ and A_l are the wavelengths parallel and perpendicular to the mean field and v is 
the mean velocity fiuctuation at the scale of the corresponding component. For an energy- 
conserving cascade, v oc A_l^^^, so equation (1) implies 



The corresponding velocity scaling law along the magnetic field is oc Ay . A critically 
balanced Alfvenic cascade leads to a scale-dependent anisotropy of the velocity field. This 
anisotropy has been demonstrated with computational simulations for both incompressible 
and compressible MHD turbulence (e.g. Maron & Goldreich 2001; Cho, Lazarian, & Vishniac 
2002; Vestuto, Ostriker, & Stone 2003). 

Can MHD induced velocity anisotropy, as predicted by GS95, be measured in interstellar 
clouds? Watson etal (2004) and Wiebe & Watson (2007) have attributed the polarization 
properties of both OH masers and thermal molecular line emission to directionally dependent 
optical depths induced by MHD turbulence. More panoramic observational views of the gas 
dynamics rely on spectroscopic imaging data of atomic or molecular line emission, most 
notably, the HI 21cm fine, and the low rotational transitions of ^^GO and its isotopomers. 



X\l/VA ~ X±/v 



(1) 



A„ oc \^'/' 



(2) 
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^^CO , and C^^O. In principle, the spatial variation of line shapes and velocity displacements 
offer a proxy view of the prevailing cloud dynamics. Recovering the form of the velocity power 
spectrum or its equivalent structure function from the spectroscopic data cubes, T(x,y,v), 
is challenging, owing to the complex integration of the velocity and density fields along the 
line of sight and the effects of line excitation and opacity that may filter or mask dynamical 
information from some fraction of the volume (Brunt & Mac Low 2004; Ossenkopf etal 2006). 

Principal Component Analysis (PCA) is a powerful method to examine spectroscopic 
imaging data of interstellar clouds. It reorders the data onto a set of eigenfunctions and 
eigcnimagcs (Hcyer & Schloerb 1997; Brunt & Hcycr 2002). Characteristic velocity dif- 
ferences, 5v, and spatial scales, r, are derived for each principal component for which the 
signal variance is distinguished from the statistical noise of the data. The set 5v,t points 
can be empirically linked to the true velocity structure function using model velocity and 
density fields (Brunt & Heyer 2002; Brunt etal 2003). The method has been appfied to a 
large set of ^^CO and ^^CO imaging observations of giant molecular clouds located within 
4 kpc of the Sun to establish the universality of turbulence within the molecular interstel- 
lar medium (Brunt 2003; Heyer & Brunt 2004). However, these studies did not consider 
velocity anisotropy. The eigenvectors were derived from the covariance matrix that was ac- 
cumulated from all spectra within the data cube with no orientation constraints. Therefore, 
any dynamical signature of anisotropy along a given axis was necessarily diluted by isotropic 
contributions to the covariance matrix. The corresponding eigenimages identified locations 
within the projected plane where velocity differences can occur but at any angle. 

In this paper, we describe a modified application of PCA on spectroscopic imaging data 
to recover structure functions along perpendicular axes (§2). In §3, the utility of this analysis 
is demonstrated on spectroscopic data cubes derived from model velocity and density fields 
from decaying MHD simulations covering a range of magnetic field strengths. In §4, we 
apply this analysis to ^^CO J=l-0 observations of a sub- field within the Taurus Molecular 
Cloud to show that such anisotropy is present and that the degree of anisotropy provides a 
coarse estimate to the strength of the magnetic field in this region. 

2. Description of Analysis Method: Axis Constrained Covariance Matrix 

To examine the degree of velocity anisotropy in interstellar clouds, we have made a 
simple modification to the application of Principal Component Analysis. A directional con- 
straint is imposed on the eigenvectors by calculating the covariance matrix from the sequence 
of spectra along one spatial axis (position- velocity slices of the data cube). The position- 
velocity image can be extracted one slice at a time to preserve spatial resolution or can be 



-4- 



generated by averaging contiguous slices to increase the signal to noise ratio. For a given data 
cube, T{x,y,v), with dimensions nx,ny,nv, the position- velocity slice along the x direction, 
averaged over thickness A in the y direction, is 

1 

Wy{x,v)^-Y_,T{x,yj,v) (3) 

where A — j2 — jl + 1. The covariance matrix for this position-velocity shce, C^, has 
components 

Cti^-J2w{x„v,)W{x„vi), (4) 
^ i=i 

(suppressing the y subscript on W) . The eigenvalue equation is solved for this axis-constrained 
covariance matrix, 

= XxUx (5) 

to produce the set of riy eigenvectors, Ux{v), that describe velocity differences exclusively 
along this particular position-velocity slice. To spatially isolate where these differences occur 
for each component, the spectra are projected onto the corresponding eigenvector, 

h{xi) = ^ W{xi, Vk)ux{vk) (6) 
fc=i 

The eigenprojection, Ix{x), has dimensions Ux x 1- The characteristic velocity difference, 
Svx, and scale, r^, are determined from the scale length of the normalized autocorrelation 
functions of Ux{v) and Ix{x) respectively (Brunt & Heyer 2002). Typically, only 4 to 5 
{Sv,t) pairs can be extracted from the axis constrained eigenvectors and projections for a 
given position-velocity shce owing to the limited spatial dynamic range. These steps are 
repeated for all position-velocity slices {jl — 1, 1 -|- A, 1 -|- 2A, Uy — A) in the data cube 
to produce a composite set of {Svx,Tx) pairs derived from n^/A sets of eigenvectors and 
eigenprojections. Similarly, to examine structure along the y axis, these steps are applied to 
position- velocity slices along the y direction averaged over x— thickness A — i2 — il + 1, 

Wxiy,v)^-J2T{^i,y,v) (7) 

i=il 

with covariance matrix 

Cli^-Y.W{yj,v,)Wiyj,vi) (8) 
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A corresponding set of {Svy,Ty) pairs are derived from Ux/A. sets of eigenvectors, Uy{v), and 
eigenprojections, Iy{y)- For each axis, we consolidate the r vahies into one pixel wide bins 
and calculate the mean and standard deviation of 5v values for each bin. Power laws are fit 
to each set to derive a relationship between the magnitude of velocity differences in the line 
profiles and scale over which these differences occur, when constrained to each axis. 



The PCA scaling exponents, a^, Oiy, are empirically linked to the scaling exponents of the 
first order velocity structure function 



(Brunt etal 2003). The first-order structure functions, 6vx = fo,a;Tj"' and 6vy = Vo,yTy^ 
provide equivalent information to the power spectrum of the velocity field along the kx and 
ky axes, respectively, for k^ — 0. Averaging over A in y (or in x) is equivalent to integrating 
along ky (or kx). 

This method offers a tool to derive velocity structure functions along any two perpen- 
dicular (spatial) axes of a spectroscopic data cube. With apriori knowledge of the local 
magnetic field direction, one could simply rotate the data cube to align the x-axis along 
this direction and determine the parallel and perpendicular structure functions. However, 
this orientation may not necessarily correspond to the angle at which velocity anisotropy is 
largest. A more rigorous test of MHD-induced anisotropy is the demonstration that velocity 
anisotropy is maximized when one of the two orthogonal axes lies along the local magnetic 
field direction. To determine the angle of maximum anisotropy, 6 max, the spectroscopic 
data cube is rotated through a sequence of angles, 9, in the plane of the sky from which the 
X and y-axis structure functions are calculated for each angle. To compare with polarization 
observations that measure position angles east of north, we define 9 as the angle measured 
counter-clockwise from the y axis. To quantify the difference between the x and y-axis struc- 
ture functions for each angle, we consider two separate measures of anisotropy. The first 
anisotropy index, ^i, is motivated by GS95 who predict differences in the scaling exponents, 

*i = l:!i—(y (11) 

Based on the results shown in §3.1 and §3.3, the second anisotropy index, ^2, measures the 
difference between the normalization constants, Vo,x5^o,y 



< 5vx >= Vo,xrx'' 

< SVy >^ V^^yT^y . 



(9a) 
(9b) 



7 = 1.69q; - 0.54 a < 0.67 




7 = 0.93a -0.03 a > 0.67 



(10b) 



*2 



- V, 



(12) 



Vo,y + Vo ,x 
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For an isotropic velocity field, and ~ 0. 

The modulation of by position angle enables a more accurate determination of the 
amplitude and angle at which the velocity anisotropy is maximized. This modulation involves 
rotation about an axis so there is degeneracy for angles 9 and 9 + 180. We find that the 
function 

■^{9) = ■^oCos[2{9 - 9max)] (13) 

provides a reasonable fit to the variation of the anisotropy index. The coefficient gives 
the amplitude of the anisotropy and the phase, 9 max, is the angle of maximum anisotropy 
that can be compared to the local field direction, {9b)- 



3. MHD Simulations 

To demonstrate that the analysis described in §2 can indeed recover the spatial statistics 
of the velocity field, we have analyzed a set of computational simulations of decaying MHD 
turbulence from Ostriker, Stone, & Gammie (2001). The models span a range of magnetic 
field strengths parameterized by the ratio of thermal to mean-field magnetic energy densities, 
/3 = c1/v\ where Cg is the sound speed for H2 and va is the Alfven velocity based on the 
mean field, (B) /y/Anp. In these decaying-turbulence simulations, the initial velocity field 
is identical for all models, so any subsequent differences in the velocity field (including 
anisotropy) arise due to magnetic effects. For each (3 model (/? = 0.01,0.1, and 1), we 
examine two snapshots at times t in units of the sound crossing time, tg- The snapshots are 
chosen such that the kinetic energy, or sonic Mach number — '^rms/cs, is comparable 
for all /3. A summary of the simulation snapshots is listed in Table 1. The cloud models 
are initially threaded with a spatially uniform magnetic field along axis 1 of the volume 
(Bo = (-Bi,0,0)); owing to periodic boundary conditions adopted for the simulations, the 
mean field (B) = Bo at all times, although the total magnetic field strength |B| changes in 
time. Depending on the strength of the mean magnetic field. Bo-, turbulent fiows can distort 
the magnetic field lines. A coarse estimate to the large scale magnetic field alignment is 
provided by the variance of each component Si, B2 and normalized by the total field 
B = {Bl + Bl + Biy/'^. For the strong field simulation {(5 = 0.01; snapshots B2,B3), 
the variance of the magnetic field components on each axis is small (< 3%), indicative 
of a spatially rigid field well ahgned along axis 1 for most cells within the volume. For the 
intermediate {(3 = 0.1; snapshots C2, C3) and weak field (/3 = 1; snapshots D2,D3) cases, the 
fiuctuations of the field values are significantly larger (10-40%), reflecting localized tangling 
of the magnetic fleld. 
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3.1. Direct Measurement of Velocity Anisotropy 

The MHD velocity anisotropy is manifest by the spectral properties of the velocity field 
along axes parallel and perpendicular to the magnetic field. Cho, Lazarian, & Vishniac (2002) 
and Vestuto, Ostriker, & Stone (2003) examined the spectral slopes of directional power 
spectra or equivalently, the 2nd order structure function, of velocity fields from computational 
simulations. Both studies found steeper spectral slopes and smaller normalization constants 
for structure functions extracted along the magnetic field direction relative to those along an 
axis perpendicular to the field. That is, the velocity field contains more power when k± ^ k 
and k\\ ^ than when k\\ ~ k and k± ~ 0, for a given k. 

To quantify the velocity anisotropy in the simulations used in this study and to compare 
with simulated observations shown in §3.2, the true 2nd order structure function, 5*2 (r), is 
calculated directly from each model velocity field, V2, 

^2(r|,,r^) = (h(x)-t;2(x + r)f) (14) 

where r = r||e|| + Tj_ej_; ey, e^ are unit vectors parallel and perpendicular respectively to 
the local mean magnetic field direction, and the angle brackets denote a spatial average over 
the volume (Cho, Lazarian, & Vishniac 2002). Here, we restrict our analysis to the projected 
plane appropriate for the velocity field to facilitate comparison with model observations in 
§3.2. In this case, the V2 component projects into the plane defined by axes 1 and 3. 
Figure [1] shows the 2nd order structure functions, 5*2(0, r^) and S'2(r||,0). Power laws are fit 
over the pixel range 5-15 to exclude the steep component at small scales that results from 
grid-scale numerical dissipation of the simulation. The amplitudes and spectral indices of an 
equivalent, first order structure function, (5*2)^/^, are listed in Table 2. Velocity anisotropy 
is clearly identified in the B2 and B3 simulation snapshots as the slope and amplitude of the 
orthogonal structure functions are different. For the intermediate (snapshots C2,C3) and 
weak (snapshots D2,D3) B-field cases, the structure functions are statistically equivalent, 
indicative of globally isotropic velocity fields with slopes (~0.5) that are typical of strongly 
supersonic, super-Alfvenic turbulent fiows. The absence of velocity anisotropy results from 
the local distortions of the magnetic field that dilute any signature to large scale anisotropy. 



^ To distinguish between the 3 spatial coordinate axes of the models and the observed projected axes, we 
reference the 3 spatial axes of the model fields as 1,2,3 and label the projected, observed axes as x and y. 
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3.2. Model Spectroscopic Data Cubes 

Observers do not directly recover the 3 dimensional velocity fields. Wide field spectro- 
scopic imaging measures line intensity as a function of position on the sky and velocity along 
an axis. The precise shape of a line profile is dependent on density, the projected veloc- 
ity component, temperature, and chemical abundance that are integrated along the line of 
sight and affected by line excitation and opacity. To place the model velocity and density 
fields from the computational simulations in the same domain as observations, we generate 
synthetic line profiles of ^^CO and ^^CO J=l-0 emission. Details of the line excitation and 
radiative transfer calculations are described by Brunt & Heyer (2002). The assumed abun- 
dance values of ^^CO and ^'^CO relative to H2 are 1.0x10"^ and 1.25x10^® respectively. 
We adopt a uniform kinetic temperature of 15 K, which corresponds to a one dimensional 
sound speed of 0.22 km s^^. The adopted mean volume density of H2 is n = 1000 cm~^. 

The choice of constructing synthetic profiles of the high opacity ^^CO emission is mo- 
tivated by two factors. First, the ^^CO J=l-0 line is the most common tracer of cloud 
structure so there are many observational data sets available to compare with these models. 
To be sure, ^^CO does not effectively probe the high density cores of molecular clouds where 
star formation takes place. However, these regions comprise a small fraction of the cloud 
mass and volume (Heyer, Ladd, & Carpenter 1996; Goldsmith etal 2008). Brunt & Heyer 
(2002) examined the effects of fine opacity on the gas dynamics perceived by observations. 
With the exception of micro-turbulent velocity fields, they found that ^^CO measurements 
reliably recover the velocity field statistics. Although the local optical depth can be large 
within a volume, the macro-turbulent velocity fields provide an effective large velocity gra- 
dient condition that allows most photons from the surface of the local volume to escape. In 
addition, owing to radiative trapping, ^^CO is detected over a broader area than the lower 
opacity lines so there are simply more measurements and information on the largest scales. 
Nevertheless, to re-examine the effects of line opacity, we also generate and analyze synthetic 
profiles of the ^^CO J=l-0 transition. 



3.3. Axis- Constrained PCA Applied to Model Data Cubes 

The utility of the analysis described in §2 is assessed by its application to the synthetic 
^^CO and ^^CO data cubes constructed from the MHD model density and velocity fields. 
Does the analysis recover velocity anisotropy when this is present in the raw velocity field, 
for the case of strong magnetic fields? Does the method verify isotropic velocity fields in the 
intermediate and weak field cases? 
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Here, we examine the synthetic spectroscopic data cubes derived from the f 2 velocity 
field for position angle 6=90°, i.e. corresponding to alignment of the x-axis with the mean 
magnetic field direction in the plane of the sky. The results of the axis-constrained PCA 
method (with A=2), as applied to all model snapshots, are shown in Figure [2J Magnetically 
aligned anisotropy is clearly identified for the B2 and B3 simulation snapshots as a separation 
of the set of points {< 6v >, r) derived respectively along the x and y axes of the model data 
cubes. This separation of points is qualitatively similar to the corresponding true structure 
functions calculated directly from the velocity fields that are shown in Figure [TJ For the 
intermediate- (C2,C3) and weak- (D2,D3) magnetic field snapshots, there is a strong overlap 
of points (< 6v >,t) derived for the orthogonal axes. This indicates velocity isotropy with 
respect to the mean magnetic field, and is in agreement with the true velocity structure 
functions for the models. 

To assess the method quantitatively, bisector fits of power laws with parameters, a, Vo, 
are fit to each set of points for each axis over the range 3 < r < 30 pixels. The scaling 
exponents, 7|j and 7_l, of the structure function are derived from the fitted parameters, a\\ 
and a±, according to equation 10. The results for the ^^CO and ^^CO model data cubes 
are summarized in Table 3. With the exception of the C3 model data cube, there are no 
significant differences between the power law parameters derived from ^^CO and ^^CO model 
cubes, demonstrating that opacity effects do not significantly skew the derived velocity field 
statistics. 

For the strong field simulations, the separation of points in Figure [2] is due to a combi- 
nation of a larger normalization constant and shallower index for the perpendicular structure 
function. Moreover, the anisotropy is stronger in the later stage simulation (comparing B3 
with B2). There are, however, discrepancies between the values of 7 determined directly 
from the velocity field in Table 2 and those determined by PCA that are listed in Table 3. 
The root-mean-square difference between power law indices is 0.12 (23%). This discrepancy 
is due in part, to the difficulty in measuring a power law index of structure functions of 
velocity fields produced by the computational simulations that have limited inertial range 
(Vestuto, Ostriker, & Stone 2003). In addition, the PCA eigenprojection along a single axis 
tends to limit the dynamic range of spatial scales over which the power laws parameters 
are derived. Despite this discrepancy of the scaling exponents, these results demonstrate 
the ability of the axis constrained PCA eigenfunctions to show a clear signature of velocity 
anisotropy induced by MHD turbulence. 
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4. The Taurus Molecular Cloud 

The Taurus Molecular Cloud provides a valuable platform to investigate interstellar gas 
dynamics and the star formation process, owing to its proximity (140 pc) and the wealth of 
complementary data. Narayanan etal (2008) present new wide-field imaging observations of 
^^CO and ^^CO J=l-0 emission from the central 100 deg^ of the Taurus cloud complex, ob- 
tained with the FCRAO 14m telescope. The images identify a low column density substrate 
of gas that contain subtle streaks of elevated ^^CO emission aligned along the local magnetic 
field direction as determined from stellar polarization measurements (Heiles 2000). Images 
of ^^CO J=l-0 integrated intensity and centroid velocity with measured polarization vectors 
from this subfield are shown in Figure [31 These show a connection between the density 
and velocity fields. While the origin of these streaks is unknown, their rigorous alignment 
with the polarization vectors strongly suggests that the interstellar magnetic field plays a 
prominent role in the gas dynamics of this low density material. 

To assess the degree of velocity anisotropy within this sub-region of the Taurus molecular 
cloud, we have applied the axis constrained PCA method to the ^^CO data from this imaging 
survey. The precise field is described by the solid box in Figure 3. We do not consider the 
^^CO J=l-0 data since the signal is weak from this low column density sector of the cloud. 
The mean, local polarization angle, derived from 16 measurements within the field is 52°±10°. 
Assuming the polarization is induced by selective absorption of background starlight by 
magnetically aligned, elongated dust grains, this angle corresponds to the local magnetic 
field direction (Purcell 1979; Draine 2003). Figure H] shows the variation of the anisotropy 
indices, and ^^2; with position angle (measured east of north) for ^^CO data within this 
subfield of the Taurus cloud. For which considers the differences in scaling exponents, the 
fitted parameters are \l/o=0.49±0.03 and 6'max =41°±2°. For \l/2, which measures anisotropy 
based on the differences of the normalization constants, \l/o=0.56±0.03 and ^max=46°±2°. 
The angle of maximum anisotropy is within 6-11° of the local magnetic field direction and 
the mean position angle of the emission streaks of ^^CO emission. The x and y-axis structure 
functions derived at 6max=4:Q° are shown in Figure O These distributions show the same 
pattern of offsets between the parallel and perpendicular structure functions measured in 
the strong field simulation snapshots (B2,B3) shown in Figure |2j For the Taurus field, the 
power law index of the structure function derived from ^^CO along the x-axis (i.e. the 
direction aligned with the polarization) is steeper (0.81±0.05) than the index of the |/-axis 
structure function (0.34±0.06). The steeper power law along the x-axis is indicative of a 
velocity field more dominated by large scales. Similar to the model structure functions in the 
strong magnetic field cases, the normalization of the y-axis structure function, Vo,y is 0.08 
km s~^ and larger than the value of the x-axis structure function (fo,x=0.02 km s^^). Thus, 
the smooth variation of density along the presumed magnetic field is mirrored by a smooth 
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variation in the velocity, and the stronger variation in density in the perpendicular direction 
(streakiness) is mirrored by a stronger variation in the velocity. Indeed, preliminary analysis 
shows that in the direction perpendicular to the projected magnetic field, displacements 
between the peaks in integrated intensity and velocity centroids are similar with typical 
values 0.2 to 0.4 pc. 

The results shown in Figures 3,4,5 are suggestive of velocity anisotropy induced by strong 
MHD turbulence, as described by GS95 and verified by computational simulations (Cho, 
Lazarian, & Vishniac 2002; Vestuto, Ostrikcr, & Stone 2003). We note that the observed 
spectral slope parallel to the field, 7|| , is steeper than the value predicted for incompressible 
MHD turbulence by GS95 but is similar to values derived for the strong field (B2, B3) 
simulations. Velocity anisotropy could be produced by processes other than MHD turbulence. 
A systematic flow of material that is "channeled" by the magnetic field may also generate 
differences in the parallel and perpendicular structure functions. Such large scale gradients 
would produce steep spectral indices (7 > 1). However, the observed high frequency variation 
of velocities perpendicular to the field are not characteristic of such large scale shear fiows. 
Regardless of its origin, the near alignment of the velocity anisotropy with the local magnetic 
field direction demonstrates the importance of the interstellar magnetic field on the gas 
dynamics within this low density component of the Taurus molecular cloud. 



4.1. The Magnetic Field Strength in the Taurus Cloud Envelope 

Since anisotropy is only evident in models with strong magnetic fields, the identification 
of such anisotropy within observational data offers a proxy measure of the magnetic field and 
its effect upon the neutral gas (Vesuto, Ostriker, & Stone 2003). Specifically, the amplitude 
of the mean magnetic field. Bo = \ (B) | = Cg ^/Attp/P, may be estimated from values of P that 
are constrained by the observations. The measured degree of velocity anisotropy is sensitive 
to the projected component of the mean field in the plane of the sky. It is improbable that 
the magnetic field threading an interstellar cloud is aligned in the sky plane. Therefore, 
measures of velocity anisotropy provide a lower limit to the value of Bo. 

Based on our analyses, the anisotropy measured in the Taurus subfield is not as large 
as in the strong field snapshots (B2,B3), but it is larger than the anisotropy hmits for the 
intermediate-field strength model snapshots (C2,C3). Given this bracketing, we can assign 
an approximate value of [3=0.03 to the Taurus subfield as a logarithmic midpoint between 
the intermediate and strong field models. Since the observed region is within the low column 
density regime of the Taurus Cloud, we set the kinetic temperature to be 15 K and the 
mean density to be 250 cm~^. These values for the temperature and density are reasonably 
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constrained by non-LTE excitation models that match the observed ^^CO and ^^CO J=l-0 
intensities from the sub-thermally excited component of the Taurus cloud (Goldsmith etal 
2008). The magnetic field strength corresponding to these values of (3, kinetic temperature, 
and gas density is 14 fiG. As noted above, this is a lower limit on the total magnetic field 
strength since the velocity anisotropy is not sensitive to the line-of-sight component of the 
magnetic field. 

Zeeman measurements of the OH line emission from the L1544 dark cloud, located ~4 
degrees to the south-west of the subfield in Taurus, identify a line of sight field strength of 11 
/iG (Crutchcr & Troland 2000). While this value is comparable to our coarse estimate of the 
field, these OH Zeeman observations are toward higher column density material {N{H2) ~ 
2 X 10^^) than is likely present in the Taurus subfield. If this higher column density reflects 
a larger volume density and if the magnetic fleld is correspondingly compressed, the field in 
the diffuse parts of the Taurus cloud may be smaller. 

The Chandrasekhar & Fermi (1953) method offers an additional measure of the magnetic 
field strength in interstellar clouds. It attributes deviations of the local magnetic field from 
the mean field direction to linear- amplitude transverse MHD waves such that 

(SB/Bp) = \Sv\/vA (15) 

where Bp is the projection of the mean magnetic field on the plane of the sky, SB and 
\5v\ are components of the magnetic and velocity perturbations transverse to Bp, and va 
is the Alfven velocity. Assuming polarization vectors accurately track the local magnetic 
field direction and transverse velocity perturbations in the two directions perpendicular to 
B are comparable, the Chandrasekhar-Fermi method is rewritten in terms of observational 
measures, 

apoi = /(47rpo)'/V,/Sp (16) 

where Cpoi is the dispersion of polarization angles measured in radians, (T„ is the line of 
sight velocity dispersion, po is the mean density of the gas, and the factor, f, accounts for 
density inhomogeneity and line of sight integration. Ostrikcr, Stone, & Gammic (2001) and 
Padoan etal (2001) determine / paO.4-0.5 from computational simulations. The dispersion 
of measured optical polarization angles within the target field is 0.17 radians. The fine of 
sight velocity dispersion determined from the ^^CO data is 0.38 km s~^. Assuming a mean 
density of 250 cm~^ and f=0.5, the derived mean field strength is 14 //G. Thus, our PCA- 
based estimate of the magnetic field strength in Taurus also compares favorably to the value 
derived by the Chandrasekhar-Fermi method. 
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4.2. The Magnetic Support of the Taurus Cloud Envelope 

The degree to which the magnetic field can support a volume against self-gravitational 
collapse is parameterized by the mass to fiux ratio with respect to the critical value, (M/ $)crii = 
O.ie/G^/^ (Nakano & Nakamura 1978). The magnetic critical index, /x, is the ratio of the 
mass to flux ratio of a volume to this critical value, 

II = (M/$)/(M/$),Ht = 7.6 X 1Q-^^N{H2)/B (17) 

where N{H2) is the gas column density in cm"^ along field lines and B is the magnetic 
field strength expressed in yuG. Owing to projections of the magnetic field and the mass 
distribution along field lines, the observed index, fiobs, overestimates the true magnetic index. 
Assuming random orientations of the magnetic field and fiattened gas distribution with 
respect to the observer, one can derive a statistical correction to the observed value, < fi >= 
fJ^obs/'^ to assess whether a volume is super-critical {< /i > greater than 1) or sub-critical 
(< /X > less than 1) (Heiles & Crutcher 2005). 

Goldsmith etal (2008) derive the distribution of molecular hydrogen over 100 deg^ of the 
Taurus Molecular Cloud using the ^^CO and ^^CO J=l-0 data of Narayanan etal (2008). 
From the Goldsmith etal (2008) image, the mean column density within the Taurus subfield 
analyzed in this study is 1.5x10^^ cm~^. For a magnetic field with strength 14 /xG, this 
column density corresponds to an observed magnetic index of /lobs — 0.81. Applying the 
statistical correction for projections, < fx >— 0.27. This low column density subfield within 
the Taurus cloud is magnetically sub-critical indicative of a magnetically supported cloud 
envelope. Such sub-critical, low column density envelopes are expected given the exposure 
to the ambient UV radiation field that maintains a sufficient degree of ionization to couple 
the neutral material to ions. The ambipolar diffusion time scale is long with respect to 
the dynamical time of the envelope. While star formation within the high density cores 
and filaments of the Taurus cloud attest to the gravitational collapse and lack of magnetic 
support within localized regions, these occupy a small fraction of the mass and area of the 
cloud. Goldsmith etal (2008) report that 50% of the mass and 75% of the area of Taurus 
have molecular column densities less than 2x10^^ cm~^. If this column density regime is 
similar to the subfield analyzed in this study, then the Taurus molecular cloud envelope 
remains magnetically supported. 

5. Summary 

We have developed an analysis method to assess velocity anisotropy within interstellar 
molecular clouds from spectroscopic imaging observations. Such anisotropy is predicted from 



-14- 



theory of strong MHD turbulence (GS95). The utility of our method is demonstrated using 
MHD simulations with varying magnetic field strengths. Velocity anisotropy is recovered in 
models with strong magnetic fields {(3 = 0.01) oriented perpendicular to the line-of-sight. 
No anisotropy is measured in simulations with the magnetic field pressure more comparable 
to the local thermal pressure, or a few times larger {f3 = 1.0,0.1). The analysis is applied 
to ^^CO J=l-0 emission from a low density sub-region within the Taurus molecular cloud. 
We detect velocity anisotropy that is aligned within ~10 degrees of the local magnetic 
field direction. This coincidence of the field direction with measured anisotropy in small- 
scale velocity variations demonstrates a strong coupling of the interstellar field with the 
neutral gas that may result from MHD turbulent flows. Our estimate of the plane-of-sky 
magnetic field strength based on our velocity anisotropy analysis is in agreement with the 
value derived using the Chandrasckhar- Fermi method. Based on our estimated magnetic field 
strength combined with column density measurements, we find that the low-density envelope 
of Taurus, which comprises the bulk of the cloud's mass, is magnetically subcritical. 

This work was supported by NSF grant AST 0540852 to the Five College Radio As- 
tronomy Observatory. ECO is supported by NSF grant AST0507315. 
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Table 1: MHD Simulation Parameters 



Model 


/3 


t/ts 


Ms 


B2 


0.01 


0.07 


7.4 


C2 


0.10 


0.04 


7.6 


D2 


1.00 


0.05 


7.2 


B3 


0.01 


0.19 


4.9 


C3 


0.10 


0.09 


4.9 


D3 


1.00 


0.09 


4.9 
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Table 2: Structure Function Parameters: (r))^/^ = VoT'^ 



Model 7|| 7^ Vo,\\ fo,± 

"^2 0.68 0.35 0.47 1.76 

C2 0.53 0.49 0.80 0.89 

D2 0.63 0.61 0.60 0.63 

B3 0.82 0.29 0.14 1.16 

C3 0.49 0.44 0.74 0.92 

D3 0.56 0.52 0.54 0.57 
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Table 3: Structure Function Parameters Derived from Model Data Cubes 



Model 




12 


CO 






i^CO 


l'o,± 


B2 


0.46 


0.26 


0.11 


0.23 


0.55 


0.32 


0.08 


0.15 


C2 


0.67 


0.53 


0.09 


0.12 


0.57 


0.44 


0.10 


0.12 


D2 


0.62 


0.53 


0.10 


0.11 


0.53 


0.59 


0.10 


0.10 


B3 


0.61 


0.15 


0.04 


0.17 


0.62 


0.23 


0.04 


0.13 


C3 


0.38 


0.39 


0.12 


0.12 


0.34 


0.53 


0.10 


0.08 


D3 


0.49 


0.31 


0.09 


0.11 


0.61 


0.41 


0.07 


0.09 
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Fig. 1. — The second order velocity structure function, 5'2(r), along axes parallel (red tri- 
angles) and perpendicular (blue circles) to the mean magnetic field direction for velocity 
fields from turbulent simulations. Velocity anisotropy is evident in the strong field model 
cases (B2, B3) as a larger scahng amplitude and shallower index for the structure function 
perpendicular to the mean field direction. 
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Fig. 2. — The axis-constrained PCA <5v >,t relationships derived from synthetic spectro- 
scopic data cubes of ^^CO J=l-0 emission for the x axis ( along the mean magnetic field; 
red triangles) and y axis ( perpendicular to the mean magnetic field; blue circles). The 
error bars reflect the standard deviation of values within each 1 pixel wide bin of r. The 
method recovers the anisotropy intrinsic to the B2 and B3 model velocity fields and verifies 
the isotropic velocity fields of the intermediate (C2,C3) and weak field (D2,D3) models. 
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Fig. 3. — (left) Image of ■'^^CO J=l-0 emission of a sub-field within the Taurus molecular 
cloud integrated over the velocity interval 5.5 to 7.5 km s~^ and (right) image of ^^CO 
velocity centroid (Narayanan etal 2008), with overlay of optical polarization vectors from 
the compilation by Heiles (2000). The molecular line emission and velocities exhibit streaks 
that are aligned along the local magnetic field direction. The solid line box outlines the area 
upon which the axis constrained PCA method is applied. The dotted line box shows the 
area within which the polarization angles are averaged to esimate the mean magnetic field 
direction. 
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Fig. 4. — The variation of the anisotropy indices, (left) and ^2 (right), with position angle, 
9. For each index, the solid line shows the fit of equation 13 to the set of points. The dashed 
vertical lines show {6b) ± Ic" inferred from optical polarization measurements of background 
stars within the subfield. The angle of maximum anisotropy is nearly aligned with the local 
magnetic field direction, which suggests a relationship between velocity anisotropy and the 
interstellar magnetic field induced by strong MHD turbulence. 
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Fig. 5. — The PCA derived 5f , r relationship derived from the ^^CO spectroscopic data cube 
of the Taurus sub-field rotated to align the x-axis with the angle of maximum anisotropy, 
(^MAx = 46°. The red triangles are points derived along the rotated x-axis and the blue circles 
are points derived along the rotated y-axis. The error bars reflect the standard deviation of 
values in each bin. The pattern is similar to that found from the B2, B3 simulation snapshots 
in Figure [2] and suggests an important role of the magnetic field on the local gas dynamics. 



